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Abstract. In this article I briefly describe how deep radio surveys may provide a 
means to identify variations in the upper end of the initial mass function (IMF) in star- 
forming galaxies at high redshifts (i.e., z >3). At such high redshifts, I argue that deep 
radio continuum observations at frequencies >10 GHz using next generation facilities 
(e.g., EVLA, MeerKAT, SKA/NAA) will likely provide the most accurate measure- 
ments for the ionizing photon rates (star formation rates; SFRs) of normal galaxies 
since their non-thermal emission should be highly suppressed due to the increased in- 
verse Compton (IC) losses from the cosmic microwave background (CMB), leaving 
only thermal (free-free) emission detectable. Thus, a careful analysis of such obser- 
vations in combination with future ALMA and JWST data, measuring the rest-frame 
far-infrared and UV emission from the same population of galaxies, may yield the best 
means to search for variability in the stellar IMF at such epochs. 



1. Introduction 

While star-forming galaxies are currently thought to b e responsible for completely 



reionizing the intergalactic medium (IGM) by z ~ 6 (e.g. iBecker et al.ll2001l : iFan et al. 



2006), the ionizi ng flux arising from star formation in Lyman break galaxies at similar 



redshifts dBouwens et alJl2007L 12008) appears to fall a factor of >6 below the minimum 



value required to maintain an ionized IGM for a given clumping factor and escape frac 



tion under the assumption of a Salpeter stellar initial mass function (IMF). ICharvl (I2008D 
has shown that this discrepancy can be reconciled by flattening the stellar IMF to have 

a slope of 1.7 if reionization occurred at z = 9 thus increasing the ionizing photon 

rate. While the idea of a top heavy IMF at these epochs does not seem completely in- 
appropriate given that low metallicity environments will favor the production of more 
high-mass stars, trying to identify such a ariation remains difficult. Here, I describe a 
relatively simple test using existing data of star-forming regions in the nearby galaxy 
NGC6946. 



2. Data and SFR Calibrations 



A multi-wavelength star formation rate (SFR) comparison is presented for 10 star- 
forming regions in t he nearby (d = 6.8 Mpc) galaxy NGC 6949, including its mildly 



starbursting nucleus (IBall et alj|1985l) . These regions were selected due their existing 



and forthcoming mid- and far-infrared spect roscopic data collect ed as part of the Spitzer 
Infrared Nearby Galaxies Survey (SINGS : iKennicutt et al.1 120031) and the project Key 
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Figure 1. A model radio to far-infrared galaxy spectrum (see lMurphyl20 09 ) . The 
individual contributions from non-thermal synchrotron, free-free, and thermal dust 
emission are indicated by dot-dashed, triple dot-dashed, and dashed lines, respec- 
tively. 



Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel (KINGFISH; PI. R. 
K ennicutt). For a much more in depth description of this analysis the reader is referred 
to lMurphy et all ((2010) and E.J. Murphy et al. (2010, in preparation). 



2.1. Radio Data 

Radio imaging at 1.4 GHz (14" x 12'/5 b eam) comes from the Westerbork Synthesis 



Radio Telescope (WSRT)-SINGS survey (IBraun e t al. 2007). Imaging at 1.5, 1.7, 4.9, 
and 8.5 GHz (15" x 15" beam) all come from |Becy 42QQ2b- The 4.9 and 8.5 GHz 



radio data included single-dish measurements for short spacing corrections. Obser- 
vations in the Ka-band (26 - 40) were taken using the Caltech Continuum Backend 
(CCB) on the GBT which simultaneously measures the entire Ka bandwidth over 4 
equally spaced frequency channels. Reference beams are measured by nodding 1'3 
away from the source; details a bout the reference beam locations can ben found in 
Figure 1 of iMurphv etall d2010h . The average FWHM of the GBT beam in the Ka- 
band was «25"among our sets of observations, which projects to a physical scale of 
« 0.8 kpc at the distance of NGC 6946. A detailed description on the performance 
of the CCB receive r, the data reduction pipeline, and error estimates, can be found in 
Mason et all (120091) . 
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2.2. IR, UV, and Optical Data 

The 8, 24, and 850 fim data used here were included in the SINGS fifth data release. 
Spitzer 70 fim and SCUBA 450 fim data, also included in the SINGS fifth data release, 
are not used in lieu of new Herschel data taken as part of KINGFISH. Only KINGFISH 
data which have a resolution better than the 25" beam of our GBT Ka-band data are 
used, which includes PACS 70, 100, and 160 fim imaging, along with SPIR E 250 and 



350 rn n imaging. The infrared photometry are fit by the SED models of Dale & H elou 
(|2002|) to derive total infrared (IR: 8 - 1000 fim) luminosities from each region. GALEX 
far-UV (FUV; 1528 A) and near-UV (NUV; 2271 A) im aging were taken from the 



GALEX archive. Har imaging of NGC 6946 was taken from Ferguson et al.l (1998). 



3. Calibration of SFR Diagnostics: Physical Assumptions 



To ensure that each SFR diagnostic can be compared fairly (i.e., borne from a com- 
mon IMF), c orresponding calibra tions are computed using St arburst99 stell ar popula- 
tion models ( Leitherer et al.|[l999 ). I choose a Kroupa IMF (iKroupalboOl ). having a 
slope of -1.3 for stellar masses between 0.1 - 0.5 Mq and -2.3 for stellar masses rang- 
ing between 0.5 - 100 Mq. A supernova cut-off mass of 8 Mq is also assumed. This 
exercise is then repeated by flattening the slope of the upper mass function from -2.3 to 
-1.5. For both IMF cases, a solar metallicity and continuous star formation (i.e., a fixed 
SFR = IMq yr" 1 ) are assumed. The resultant Starburst99 output relevant for calibrat- 
ing SFR diagnostics described below for both IMFs at an age of 100 Myr are given in 
Tabled] 



Table 1 . Starburst99 Output for Kroupa and Top Heavy IMFs 



IMF 


G(H°) 


qSNR 










(photons s 1 ) 


(SN century -1 ) 


(ergs s _1 Hz -1 ) 


(ergs s" 1 ) 


(ergs sr 1 ) 


Kroupa (-1.3,-2.3) 


1.37 x 10 53 


1.16 


1.14 x 10 28 


1.86 x 10 41 


2.58 x 10 4 3 


Top Heavy (-1.3,-1.5) 


7.78 x 10 53 


2.69 


3.57 x 10 28 


1.06 x 10 42 


8.36 x 10 4 3 


Top Heavy /Kroupa 


5.58 


2.32 


3.13 


5.69 


3.24 



Above Starburst99 output is taken after 100 Myr of continuous star formation 
(1 Mq yr -1 ) assuming solar metallicity. 



3.1. Radio SFRs 

Radio continuum emission from galaxies is generally composed of two optically thin 
components: non-thermal synchrotron emission associated with cosmic -ray (CR) elec- 
trons accelerated in a galaxy's magnetic field, and thermal bremsstrahlung (free-free) 
emission around massive star formation regions. The origin of each of these com- 
ponents lies in the process of massive star formation, which is most likely the phys- 
ical underpinning of the nearly ubiquitous far-infrared (FIR)-radio c orrelation (e.g. 
iHelou et al.lfl985l : ICondonlll992l : lYun et al.l l2001: Mu rphy et al. 2006a, 2008). In this 



simplistic model, the non-thermal radio luminosity of a galaxy at frequency v is re- 
lated to the SFR by supernova (SN) rate, #snr, while the free-free radio luminosity is 
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directly proportional to the number of ionizing photons, <2(H°), along with a kinetic 
temperature (T e ) dependence, such that, 

SFR V _ NT oc ^snr oc v a L^ T ; SFR V _ T oc Q(H°) oc v 01 Tf A5 L T v , (1) 

respectively, where a is the non-thermal spectral index. Thus, the SFR derived from 
the total radio emission at any one frequency can be described as a combination of each 
term such that, 

SFR v ocL v /(v- a + y- 01 r e - 45 ). (2) 
3.2. IR, UV, Ha, and Hybrid SFRs 

Similarly, SFRs can be expressed using the output of Starburst99 in terms of IR, UV, 
and optical (Ha) observations, in a much more straightforward manner. IR, UV, and 
Ha SFR estimates go as, 

SFRircxLfr; SFRuvcxLijv; SFR H( * oc Q(H°) oc L Ha . (3) 

Since UV and optical wavelength suffer significantly from extinction, making it diffi- 
cult to properly convert such measurements into an accurate SFR, a number of hybrid 
methods have been introduced into the literature. The simplest is using a combina- 
tion of IR and (observed) UV SFR estimates to account for both the obscured and 
unobscured components, resulting in a 'bolometric' or 'total' SFR. Another hybrid di- 
agnostic, which probably has a more solid physical backing, is the linear combination 
of observed Ha and 24 //m observations t o obtain an extinction corrected SFR (e. 



oi observed tla ma z4 //m obs ervations to obtain an extinction corrected af K (e.g ., 
Kennicuttetalll2007l : ICalzetti et al.ll2007l : iKennicutt et aDl2009l : ICalzetti et al.l l2010): 



this diagnostic has been empirically calibrated against extinction corrected NIR re- 
combination lines. Accordingly, these hybrid SFRs have the following observational 
dependencies: 

SFRtot x L m + L uv ; SFRmix oc g(H°) oc L Ha + b x L 24 m , (4) 



where b is an empirically derived coefficient (for a detailed discussion see ICalzetti et al.l 
2010h . 



4. Comparison of SFR Diagnostics 

In Figure |2] SFRs from each region in NGC 6946 are plotted using the various calibra- 
tions discussed in §3: Ha+ 24 jum SFRs are plotted against the thermal, non-thermal, 
and total 33 GHz SFR calibrations in the top row of panels, while the IR+UV SFRs 
are plotted against the same 33 GHz SFR calibrations in the bottom row of panels. I 
assume that the thermal radio continuum at 33 GHz provides the most accurate SFR 
estimate among all diagnostics as it is the most straightforward to relate to the ionizing 
photon rate (assuming that the spectral decomposition into thermal and non-thermal 
components is done correctly). SFR estimates using the thermal radio continuum at 
33 GHz and the combination of Hor+24 //m emission appear to agree quite well over 
then entire luminosity range, having a median near unity, albeit, with some scatter. 
The non-thermal 33 GHz SFRs appear to underestimate the actual SFR by a factor of 
~2, on average, while showing good agreement in the nucleus. Also notable is the in- 
creased scatter among the extranuclear Hn regions. This result can be attributed to the 
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Figure 2. A comparison of SFR diagnostics discussed in §3. From left to right, the 
top row shows the differences between the Ha+24 fim and the 33 GHz thermal, non- 
thermal, and total SFR estimates. The same is shown in the bottom row, except the 
33 GHz SFRs are compared to the IR+UV SFRs. In the top left corner of each panel, 
the median of the ordinate/abcissa ratio values, along with the dispersion about the 
fit to each trend are given. 
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Figure 3. Same as Figure 12 except SFR calibrations have been derived using a 
top-heavy IMF (see Table[T|). 
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fact that the population of CR electrons in the extranuclear star-forming regions may 
be typically older, on average, than those being freshly injected into the ISM near the 
starbusting nucleus. Thus, they have had time to propagate much larger distances, and 
no longer contribute to the non-ther mal emission in the immediate vicinity of the Hn 
regions dMurphy et alj|2006bl . l2008h . On the other hand, the SFR estimated using the 
total radio luminosity (i.e., the combination of the free-free and non-thermal emission 
as given in Equation 2) shows good agreement with the Ha+24 /im estimate. This is 
due large thermal fractions at 33 GHz. The IR+UV SFR estimate also appears to do 
generally well relative to the 33 GHz thermal radio continuum, save the nucleus, which 
has IR+UV SFR which is ~2.5 times larger. 

In Figure [3] the same comparison of SFR diagnostics is shown as in Figure |2j ex- 
cept that the top heave IMF calibrations have been used (see Table [T]). The Starburst99 
output suggests that the ionizing flux will increase by a factor of «6 by moving to a 
top heavy IMF, while the UV and IR output will only increase by a factor of «3. One 
of course expects that SFRs which depend only on the ionizing flux should be similar 
independent of the IMF. However, based on these differences in the calibrations, SFRs 
proportional to the ionizing flux should be a factor of ~2 larger those which depend on 
the UV and/or IR luminosities for a top heavy IMF compared to a standard Kroupa IMF. 
In looking at the top and bottom left panels, this is observed; the thermal radio contin- 
uum and Ha+24 /urn SFRs diagnostics still yield a nearly one-to-one correlation with 
one another while the IR+UV SFRs yield a SFR which is a factor of ~2 discrepant from 
the 33 GHz thermal radio continuum estimate (i.e., a factor of ~2 larger since the IMF 
of NGC 6946 more like a Kroupa IMF). Similarly, the same factor of ~2 discrepancy 
is found when comparing the IR+UV SFRs to the SFRs estimated from the total radio 
continuum. I note that all of these three quantities used to estimate the SFRs in this 
case are as observed; no extinction correction was necessary, nor does it depend on an 
empirical calibration. Thus, by having UV, IR, and thermal radio continuum emission 
in hand, one should be able to quantify such discrepancies between SFR diagnostics, 
and search for sources with potentially top heavy IMFs. 



5. Radio Emission from Galaxies at High-z 



For a typical star-forming galaxy, CR electrons prima rily lose their energy due to syn- 
chrotron and inverse Compton (IC) proces ses (e.g.. ICondonlll992h . although escape 
also plays a role (e.g jHelou & Bicay|[T993l) . At z = 0, U C mb ~ 4.2 x 10" 13 erg cm -3 , 



which is significantly smaller than the radiation field energy density of the Milky Way 
(i.e., Umw ~ 10~ 12 erg cm -3 ). Thus, CR electron energy losses from IC scattering 
off the CMB are negligible at low redshifts. However, Uqmb °c (1 + z) 4 , making 
such losses increasingly important with redshift. As an example, by z ~ 3, Uqmb ~ 
1.1 x 10~ 10 erg cm -3 ; equating this to the magnetic field energy density Ub = B 2 /(Sn) 
results in a corresponding magnetic field strength of ~50fiG. This is nearly an order 
of magnitude larger than the ambient field strength in the solar neighborhood. Con- 
sequently, the non-thermal component of a galaxy's radio continuum emission will be 
increasingly suppressed with increasing redshift, eventually resulting in only the ther- 
mal component being detectable. 

It is worth noting that additional energy-loss terms may become increasingly im- 
portant in the case of galaxies hosting strong starbursts. In such systems, whose energet- 
ics and ISM may be vastly different, ionization, bremsstrahlung, and adiabatic cooling 
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Figure 4. The relative cooling fractions for synchrotron, bremsstrahlung, and 
ionization processes for 1.4 GHz emitting CR electrons as a function of in- 
creasing magnetic field strength under the assumptions that (v/GHz) = 1.3 x 
l(T 2 (B///G)(£/GeV) 2 , (n ISM /cm- 3 ) ~ (fi/10 //G) 2 , and U mA ~ U B (i.e., internal 
IC and synchrotron cooling timescales are similar). 



through advection out of a galaxy by galactic scale winds can all play a significant role 
in cooling CR electrons. This is illustrated in Figure 0] where I plot the fractional energy 
budget for 1 .4 GHz emitting CR electrons as a function of magnetic field strength (see 
caption for assumptions). Accordingly, CR electron cooling via processes other than 
synchrotron emission may become even more efficient if the large scale magnetic fields 
approach and exceed mG strengths. However, other physics may also be at work in 
such configurations. For example, additional synchrotron emission from secondary e ± 
arising from n and decay as CR nuclei inelastically scatter of the interstellar gas, 
may play a role in produ cing additiona l synchrotron emi ssion to compensate for these 
additional energy losses (Murphyl l2009l : lLacki et al.ll201fj|) . This will not be the case for 
normal star-forming galaxies. 



6. Conclusions 



I have investigated a simple test which may prove useful for identifying IMF variations 
in galaxies at high z (i.e., z > 3). Using a combination of high frequency (>10 GHz) 
radio data along with far-infrared and UV data, one can compare radio vs. IR+UV SFR 
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estimates which should become increasingly discrepant as the IMF shifts towards one 
that is top heavy. This test makes use of the fact that radio continuum emission becomes 
dominated by thermal (free-free) emission at higher frequencies, and in addition, that 
the non-thermal emission of high-z galaxies should become increasingly suppressed 
due to increased IC losses of CR electrons off of the CMB. Thus, deep surveys at 
frequencies such as the X-band (8-12 GHz) should detect the free-free emission from 
galaxies at z > 3 directly. Such tests to identify variations in the IMF of high-z systems 
should be possible using deep survey imaging from upcoming radio (EVLA, MeerKAT 
& SKA/NAA), mm/IR (ALMA/CCAT), and optical/NIR (JWST) facilities. 

Acknowledgments. I would like to thank the organizers of the UP2010 confer- 
ences for putting together an exciting and diverse program, as well for letting me be a 
part of it. 
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